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We propose an approach to analyze the dissipation properties of coupled cavity arrays. Employing 
a kind of quasi-boson, it is shown that the coupling to a bath renormalizes the localized mode and the 
interaction between cavities. By virtue of without having to mention the coordinates of bath, this 
approach would be great conceptual and, moreover, computation advantage. Based on the result, 
a single-photon transport in the array is examined, and the total transmission rate is presented. 
Besides, we also suggest a parameter to scale quality of the array. 
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Coupled cavity arrays(CCAs), the effective and ma- 
nipulatable many-body system are known playing 
a key role in quantum information and quantum de- 
vices progressing [2U3[, mimicking and studying strong 
correlated physics [4-6]. In the term "CCA", "cavity" 
refers to a region of space within which photon can be 
efficiently confined, and "coupled" usually implies in- 
dividual cavities resonantly coupled to each other via 
the evanescent field. The boost of experimental tech- 
niques and advances in theories, over the past years, 
have paved the way for researching CCA systems more 
in-depth. CCAs are not only explored in various forms, 
but also manufactured with smaller mode volume, higher 
quality factor, more accurate addressability, and in- 
creased number, of cavities lying in the set-ups 0, 
Moreover, the theoretical proposals, paradigmatically 
like tight-binding model(TBM) @, [nj, Bose-Hubbard 
model(BHM) @,ElHI, and effective spin model IHG3, 
have been put forwarded and predicted a large quantity 
of novel applications [LiT - tl6| . 

Despite the substantial progress, both in experiment 
and theory, still many important issues are unsolved. A 
problem recurs often and in more than one facet is how 
dissipation would behave in CCA systems @, It is 
always referred to a system, having numerous degrees of- 
freedom, interacts with a bath, with many, in principle 
infinitely many, modes. The huge Hilbert space results in 
great hurdle, even challenge, to fully describe the proper- 
ties of such system, as yet, there is still lack a method to 
understand and calculate properly. However, the surge of 
interest in, so called, quantum manipulation and quan- 
tum simulation makes the problem becoming urgent to 
be solved. 

In this paper, we address the issue and show that, for 
weakly intercavity coupled high-Q array, most experi- 
ments in CCAs are carried out such that favour con- 
ditions where coupled- mode approximation is valid [l8| . 
the problem can be cured by using a kind of quasi-boson 
picture. Which essence is discarding the coordinates of 



bath, and describing the system via a effective Hamilto- 
nian. Throughout our discussion, we base on the vac- 
uum, realized, and most general ID CCA. The reasons 
are as follows: Firstly, photon leaking from cavity mode is 
one of the main process of dissipation. Secondly, the ap- 
proach could also account other sorts of dissipation, like 
atomic decay, in CCA systems. Finally, and most im- 
portantly, it maintains the universality of the proposed 
approach. According to examine the single photon trans- 
port in CCA, the theoretical results are in well agreement 
with experiments. Which indicates that our work may 
open up wealth of possibilities for studying dissipative 
CCA systems conveniently. 

We are starting by reviewing the configuration of ideal 
CCA. As shown in Fig. 1(a), N cavities, with a single 
mode characters by frequency u> c , are arranged in a pe- 
riod L. The coupling parameter in CCA is mathemati- 
cally described by a specific overlap integral a. Because 
individual modes are confined efficiently, only modes of 
nearest neighbor cavities have a small, but non- vanishing, 
overlap, denoted by auj c . Such configuration is well un- 
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FIG. 1: Schematics of CCA. (a) Individual cavities resonantly 
coupled to each other due to the overlap of their evanes- 
cent fields of nearest neighbors. (b) The coupling of CCA to a 
bath. Each resonator has a leakage rate 7-(c) Effective treat- 
ment in quasi-boson picture,where system can be regarded as 
a chain of quasi-bosons. 
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derstood by using TB Scheme and forms the basis for 
other CCA systems pj. In the Coulomb gauge, these 
localized modes, labelled by <pj, obey the Maxwell equa- 
tion dHH: 



e(r)w* 



and a is given by 



tpj - V x (V x ipj) = 0, 



dr[e (r) - e(r)]<pj>j 



(1) 



(2) 



Where r is a given, in fact, three dimensional vector, 
eo(r) the dielectric constant of single cavity, and e(r) the 
periodic dielectric constant of array. 

When N is large enough, periodic boundary condition 
of CCA is fulfilled and the Hamiltonian (with % = 1) 
reads 



UJ, 



OLUJr 



0\i'> 



(3) 



The bosonic operator c^{cj) creates(annihilates) a ex- 
cited state at jth cavity, j/\ sums all pairs of cavities 
which are nearest neighbors. 

As a result of TB interaction, the whole system is no 
longer a monochromatic field but splits into N resonant 
modes, which are well explained as liner combination of 
individual cavity modes, and forms a narrow band in 
vicinity of uj c . The spectrum takes the form 



u)(k) = u) c + 2auj c cos k n L 



(4) 



and has been observed experimentally by measuring the 
transmission-phase properties (20| . However, it is worth 
stressing that the wave vector, k n — j^ij; for n = 1 to 
N, does not has a direct meaning in terms of photonic 
momentum but to analog of the lattice vector in solid 
state physics. 

Since we have not added any other features, like 
Jaynes-Commings interaction and Kerr interaction, the 
formalism described above is the most common founda- 
tion of CCA systems. Therefore, it is advisable to analyze 
the dissipation properties base on the scheme. In what 
follows, we will approach the problem in two steps. 

Firstly, we consider a single cavity coupled to a bath 
composed of the infinite set of harmonic oscillators [2l[ . 
The bath generally has a continuous spectrum character- 
ized by u) r and the density of states described by p(ui r ). 
For simplicity, assuming here only one excited state oc- 
cupied by either cavity or bath, and the corresponding 
probability amplitude denoted by e c and e r respectively. 
Besides, because the totality is conservative, one can 
write the eigenvalue equation as 



H\cp) = to\ip), 



(5a) 



with 



H = uj c CqCo + / dui r r + r + / duj r [r]* (uJ r )r + cq + h.c. 



(5b) 



\ip) = e c Co"|0) + / duj r p{uj r )e r r 



(5c) 



r + (r), which satisfies the commutation relation 
[r(w r ), r + (uj' r )] — S(uj r — v' r ), creates(destroys) an 
excited state of bath. rj(ui r ) represents the coupling 
strength between the two, to the total energy, and |0) 
the vacuum state. 

Taking the inner product first with Cq |0) and then with 
r+|0) to Eq.(5a), 

w c e c + / duj r p(uj r )r](uj r )e r — uie c , (6a) 



From Eq.(6b), e r 
Eq.(6a), 



uj r e r + rj* (uj r )e c — cue r . (6b) 
e c , and plugging it into 



V 0"r) , 

Ld—UJ r 



. , . .\T](uj r )\ 2 

w c e c 4- / dui r p(us r ) e c 



: uie c 



(7) 



Note that 



du! r p(u!r) 



dui r p(uj r ) 



\rj(ujr 



uj — uj r + iS 



P I dw r p(u) r ) — r ^ i7rp(w)|?7(ti;)| 2 (8) 



oj — iv r 

In the above derivation, we have extended the in- 
tegration into complex plane and used the relation 
lim — pr- = P- — inS(x), where P denotes the Cauchy 

y _ )0+ x+iy x V )i J 

principal value, x and y are real variables. 

It is reasonable to assume that uj can cause excitation 
of the cavity sharply peaks around w c . Therefore, we can 
evaluate Eq.(8) at uj = uj c , 



P I duj r p{w r ) — r ^ ~ P f dui r p(u! r ) — r ^ 



iirp(oj)\i](oj)\ 2 w iirp(oj c )\r)(uj c )\ 2 = ij. 



Suj c , 
(9a) 

(9b) 



duj c is known analogous to the Lamb shift and signif- 
icantly small in the case of coupling to surroundings 
weakly. 7 is the decay rate, which indicates a finite life- 
time of cavity mode. 
Thus Eq.(7) becomes 



(uj c + Slu c — i^)e c — coe c 



(10) 



Which means that due to the coupling, the cavity mode is 
renormalized by reckoning in frequency shift and intrin- 
sic loss. The above expression is equivalent to Eqs.(6a) 
and (6b), however, does not contain degrees of freedom 
of bath. It motivates us to introduce a quasi-boson 
described by b and having a complex eigenfrequency 



3 



uj e ff = uj c — ry, where 6lu c has been absorbed into 
lj c , to rcdescribe the cavity mode. And then rephrasing 
Eqs.(5a)-(5c), 



H eff \tp) =U) e ff\<f), 



(11) 



with the effective Hamiltonian H e ff — u e ffb + b and now 
\tp) = e c b + \%) referred to as quasinormal-mode [HI . Be- 
cause of loss energy, the system would nonconservative, 
and the corresponding operators non-Hermitian. To com- 
pare the two descriptions, the communication relation of 
b reads [b, b + ] = 1 + i^j-. Clearly, ^ in order of thus 
bosonic communication relation is approximately satis- 
fied. 

Then next, we return to the case, see Fig. 1(b), CCA 
coupled to a bath and each resonator has a leakage rate 7. 
It is verified experimentally the main sources of loss are 
individual cavities, while the additional loss caused by 
periodic structure is negligible Q. Combining the char- 
acteristic cavities are weakly coupled, such system can 
be regarded as a chain of quasi-bosons, see Fig. 1(c), and 
mapped safety onto TB scheme. According to Eq.(l), 
the associated eigenmodes, labelled by ipj, satisfy 

€ -^(uj 2 c + 7 2 )^ - V x (V x ^) = 0. (12) 

For 7 2 is Q 2 orders of magnitude smaller than uj 2 , the 
minimal loss on each lattice site does not generate notice- 
able alteration to localized modes. Which also illustrates 
that quasi-boson picture is an excellent approximation to 
the established mode. 

Consequently, the relevant overlap integral, a', is given 

by 



'= / dr[e (v) - e(tM^ 



+1 



dr[e (r) -e(r)]tpj(p jA 



a. 



(13) 



Hence, we reach the familiar Hamiltonian but take dis- 
sipation into account, 



H 



aUJ eff b 3 b l 

<J,3') 



(14) 



Yet interestingly, without having to mention the external 
degrees of freedom, the effective treatment would be of 
great conceptual and, moreover, computational advan- 
tage rather than treatment of universe [17] . One key fea- 
ture is now the loss seems owing to the nonideal bound- 
ary but not field oscillation, viz described by a constant 
but not operators. In addition, it should also to point 
out that the specific impact of renormalization to inter- 
action terms may vary from case to case, nevertheless 
all of those represented by a small quantity aj. This is 
consistent with the conditions discussed previous. 

To demonstrate the validity of our approach, we con- 
sider now the single-photon transport in the CCA. In 



the simplest possible context, we assume that a photon 
has somehow been inject into 1st cavity and propagat- 
ing to the right. The frequency, uj, of photon satisfy 
dispersion relation (4), thus photon hopping can occur 
between neighboring cavities due to the overlap of the 
light modes. So the problem we treated can be described 
by Hamiltonian (14). Furthermore, to focus on the total 
transmission rate, we can restrict us to solve the station- 
ary Schodinger equation 



H\ip) =u\ip), 



(15) 



with 



£e,&+|0), and take [H, Hi 



j 

a ik n sL 



+ ~ik n sL 



+ r 3 e~ iknSL s < j 



s> j 



= 1 



to N. Where Tj and tj denote the local transmis- 
sion amplitude and reflection amplitude of photon 
respectively. Solving Eq.(15) by using the continuous 
condition e.j- — ej+ at jth site and the constraint 
condition I 1 2 + \tj\ < 1 due to the irreversible loss of 
energy, we get 



k cos k n L — 7 



(7 + £| sin k n L\ — n cos k n L) — in\ sin k n L\ 



J2k n jL 



(£ — in) \ sinfc„L| 



(7 + £| sin k n L\ — k cos k n L) — ik\ sin k n L\ 



(16a) 
(16b) 



Above, £ = 2ato c and k = 2aj for compactness, e l2 n ^ L 
is position-dependent global phase but does not affect the 
transport properties, and the absolute value sign is need 
for energy conservation. 

Before proceeding, here we briefly outline some of the 
main features of Tj and tj. The nonzero reflection am- 
plitude is caused by local loss. Under the circumstance 
of system is confined in one dimension, incoming photon 
having possibility to escape toward the opposite direc- 
tion. Note however, this possibility would not make pho- 
ton enters the previous cavity and becomes left-moving 
photon, but eventually decay to other dimensions. Local 
loss also leads to the nonunitary transmission amplitude. 
By dropping the second-order small quantity k, the max- 
imum of transmission coefficient approximates jf^j^i, 
which means the local transport properties is determined 
by the competition between photon hopping and decay, 
since they are the only channels photon can leave a cer- 
tain cavity. 

And then, the total transmission rate, T, can be intu- 
itively written as 



T 



life 



27V 



(17) 



The transmission spectrum, shown in Fig. 2(a), retains 
the symmetry of dispersion relation (4) and vanishes at 
band edges. When the propagating photon is on resonant 
with individual cavity, uj = uj c , the spectrum exhibits the 
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FIG. 2: (a) Transmission spectrum for single-photon trans- 
port in CCA. The total transmission rate of each splitted 
resonant mode is denoted by dots and fitted by a curve. Fol- 
lowing the experimental parameters in Ref. [jj, Q = 1.1 X 10 f) 
and £ = 6.47 x 10 -4 cj c , the output power, concretely like when 
N = 60, is well agreement, (b) Dependence of the maximal 
transmission rate on quantity factor (. In Fig. (2a), £ equal 
to 71.17, 11.86, 7.117, and 3.559 respectively. 



maximum, T„ 



While the ratio between 



7 (i+7/£) 2JV ■ 

local loss rate and intercavity coupling strength is far less 
than one, we can take (1 + f) 2N = 1 + 2^2 -| , thus 



Tn 



1 



1 



(1 + N 1 /if (1 + N/aQf 



(18) 



By substituting Q = and £ = 2aw c , the maximal 
transmission rate now is described directly with three 
essential parameters of CCA. It is helpful to define a 
new quality factor, £ = to scale CCA's transport 
properties, which lead to 



T„ 



(i + C) 5 



(19) 



Furthermore, the transport loss mainly stems from 
cavity-mode decay, thus £ could reflect as well as dis- 
sipation properties for other CCA systems. A high-£ ar- 
ray, see Fig. 1(a) for example, is often referred to steep or 
sharp spectrums. 

In summary, to aim at the descriptive difficulty caused 
by the coupling of CCA systems to environment, we have 
proposed a kind of quasi-boson picture and shown its 
effectiveness by analyzing the single-photon transport. 
Here we would like to emphasize the generality of our 
approach, which is capable of treat ing dynamical prob- 
lems and other sorts of dissipation [2J, [25[ and provides 
a starting point for discussing more complicated situa- 
tions 
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